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Using the example of the elastic vrAf-amplitude, we discuss the low energy expan- 
sion of QCD amplitudes in the sector with baryon number one. We show that the 
chiral expansion of these amplitudes breaks down in certain regions of phase space 
and present a framework which leads to a coherent description throughout the low 
energy region, while keeping Lorentz and chiral invariance manifest at every stage 
of the calculation. We explain how to construct a representation of the pion nu- 
cleon scattering amplitude in terms of functions of a single variable, which is valid 
to O(q^) and properly accounts for the tttt- and TrAf-cuts required by unitarity. 



1 Introduction 

The chiral symmetry, which the Lagrangian of QCD reveals in the limit of 
vanishing quark masses, can systematically be incorporated into an effective 
Lagrangian describing the interaction of the lowest lying baryons and mesons. 
The effective description is based on a simultaneous expansion of the QCD 
Greens functions in powers of the light quark masses and the momenta of the 
Goldstone bosons. 

In the vacuum sector of the theory, where the only low energy singularities 
arise from the propagation of the Goldstone bosons, dimensional regulariza- 
tion yields homogeneous functions of the Goldstone momenta and masses, so 
that each graph has an unambiguous order in the chiral expansion. The sin- 
gularity structure in the sector with baryon number one is more complicated. 
The expansion of the corresponding amplitudes in the momenta and masses of 
the Goldstone bosons leads to an expansion of the nucleon propagator around 
the static limit. This expansion of the nucleon propagators is implemented 
ab initio in the framework called heavy baryon chiral perturbation theoryEI 
(HBCHPT). However, as we will show in section ||, it fails to converge in 
part of the low energy region. The breakdown is related to the fact that the 
expansion of the nucleon propagator ruins in some cases the singularity struc- 
ture of the amplitudes. This makes it desirable to perform the calculations 
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in a relativistic framework. In doing so, the correct analytic properties of the 
amphtudes are guaranteed and one can address the question of their chiral 
expansion in a controlled way. 

In the relativistic formulation of the effective theory a technical compli- 
cation arises from the fact that in a standard regularization prescription, like 
dimensional regularization, the low energy expansion of the loop graphs starts 
in general at the same order as the corresponding tree diagramsH. Since the 
contributions that upset the organization of the perturbation expansion stem 
from the integration region of large loop momentum of the order of the nu- 
cleon mass, they are free of infrared singularities. In section ^ we show that 
the infrared singular part of the loop integrals can be unambiguously sepa- 
rated from the remainder, whose low energy expansion to any finite order is 
a polynomial in the momenta and quark masses. Moreover the infrared sin- 
gular and regular parts of the amplitudes separately obey the Ward identities 
of chiral symmetry. This ensures that a suitable renormalization of the effec- 
tive coupling constants removes the infrared regular part altogether, so that 
we may drop the regular part of the loop integrals and redefine them as the 
infrared singular part of the integrals in dimensional regularization, a proce- 
dure referred to as infrared regularizatiora. The representation of the various 
quantities of interest obtained in this way combines the virtues of HBCHPT 
and the relativistic formulation: Both the chiral counting rules and Lorentz 
invariance and are manifest at every stage of the calculation. 

Ih the meantime this framework has been used to calculate the scalarS, 
axialB and- electro-magnetic formfactoraa as well as the elastic pion nucleon 
amplitudeQ to fourth order in the chiral expansion.^ In section ^ we discuss 
the structure and the low energy expansion of the pion nucleon amplitude. 
To fourth order in the chiral expansion the amplitude can be given in terms 
functions of a single variable, either s, t or u. 



2 Low energy singularities of one nucleon amplitudes 

The low energy singularities of the one nucleon amplitudes are due to the 
propagation of pions and nucleons. In the example of the pion nucleon ampli- 
tude these include poles at s = to^ and u = due to one nucleon exchange, 
branch points at s = (M^ -|-mjv)^ and u = (M^ -|-mjv)^ due to ttTV intermedi- 
ate states as well as a branch point aXt = 4M^ due to TTTr-intermediate states 
in the t-channel (see figure |^). Note that some of the graphs, like figs. |l|d and 
ll|e, also involve a cut for t > 4m^ due to the presence of NN intermediate 
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(a) (b) (c) (d) (e) 



Figure 1. Examples of contributions to the nN scattering amplitude with different singu- 
larity structure. 



states in the t-channel. This singularity Ues however outside the low energy 
region and is absent to any finite order of the chiral expansion. 

In performing the low energy expansion we must choose the kinematic 
variables in such a way that the singularities of the amplitude stay put in 
order to avoid a breakdown of the expansion near the singularity. We will 
illustrate this first with the trivial example of the nucleon Born term and 
then discuss the singularity structure of the triangle graph, figure ^c. In 
addition to the cut for t > ^M^, this graph has a branch point at 

t = 4Ml - % 
ml 

on the second Riemann sheet. Since the difference between the two singular- 
ities is of higher order, they will coalesce if one performs the chiral expan- 
sion, causing a breakdown of the expansion in the vicinity of the threshold at 
t = 4M2. 



2.1 Born term 

Denoting the momentum of the incoming nucleon and pion with P resp. q, the 
momentum dependence of the s-channel Born contribution to the amplitude, 
fig. ^ is given by 

^fi °^ 5 TTT"; — — ■ (1) 

ml-{P + qY 

If we expand in the momentum of the incoming pion, the propagator denom- 
inator gets replaced by the geometric series 



ml-{P + qY 2P-qy^ 2P ■ q 

The location of the singularity has moved to P • g = or s = ml + and 
the expansion breaks down for \P ■ q\ < Ml/2. Clearly, this problem can 
be cured by performing the expansion in q at fixed P = P + q/2 instead of 
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fixed P. Note however that we need to choose appropriate variables for each 
graph separately in order not to spoil its specific singularity structure. In the 
example of the graph shown in fig. ^ the expansion at fixed P is convergent, 
while the expansion at fixed P breaks down in the vicinity of the threshold 
at P ■ q = tti^Mt^. 



2.2 Triangle graph 

The discontinuity of the triangle graph, fig. |l|, can be obtained from the ex- 
tended unitarity relation 

luvfi(t) - {1 - ^Ml/ty t'j[tf flit) , 4Af2 < t < IMll . (3) 

The quantities [t) represent the i-channel partial waves of the ttN scattering 
amplitude. The quantum number J stands for the total angular momentum, 
while the lower index refers to the spin configuration . There is no need to in 
addition specify the isospin quantum number, because it is determined by the 
total angular momentum: J even ^ / = 0, J odd ^ / = 1. The functions 
t^jit) denote the partial wave projections of the tttt scattering amplitude. To 
order O(g^) the evaluation of this condition is particularly simple, because 
the TTTT scattering amplitude is needed only to leading order, where it is a 
polynomial. Hence only the partial waves with J = 0, 1 contribute. 

To obtain the imaginary part of the triangle graph, we insert the lowest 
partial waves of the nucleon Born term into the unitarity relation (^. For 
J — the pertinent partial waves are 







The function /(k) stands for 



_ arctauK _ yj[t -AM"^) (4rn^ - t) 

''^ i-2Af2 • 

The problem addressed above shows up in this expression: The quantity k 
represents a term of order 0{l/q) and the standard chiral expansion of /(k) 
corresponds to an expansion of arctauK around k = oo. This expansion 
only converges for |k| > 1. In the vicinity of the threshold at i > 4M^ this 
condition is not met and the chiral expansion diverges. The breakdown of the 
expansion is related to the fact that arctauK has a branch points at k = ±i. 
To 0{q'^), the imaginary part of the invariant amplitude Z?+ (see eq. (Q) for 
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its definition) can be obtained from 



ml 



The analytic continuation of the scattering ampUtude to the second Rie- 
niann sheet therefore contains a branch point at 

{t - 4 Ml) {Ami ^_i^t = 4Af2 - 
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3 Infrared regularization 

Using the example of the self energy integral, we explain how to separate 
the infrared singular part of the loop integrals from their regular part. As is 
evident from eqns. (|^) and the nucleon propagator counts as a quantity 
of 0(1/(7) at tree level. While the infrared regular part of the loop integrals 
violates this counting, we show in the following that the infrared singular part 
respects the counting rules established at tree level. 



3.1 Self-energy integral 

The simplest example of a one loop contribution, the self energy graph, is 
shown in figure pk. The corresponding scalar loop integral has the form 



P-k 
(a) 



' J J J 

(b) 



Figure 2. a.) Self energy graph, b.) Form of the general one loop graph. 



^ ' i J {2TrY M2 - /fc2 ^2^ _ _ ;,)2 ^ J (277)-^ a b ' 

We need to analyze the integral for external momenta in the vicinity of the 
mass shell: P — m,siV + r, where w is a time-like unit vector and r is a quantity 
of order q. In the limit M 0, the integral develops an infrared singularity, 
generated by small values of the variable of integration, k — 0{q). In that 
region, the first factor in the denominator is of 0{q'^), while the second is of 
order 0{q). Accordingly the chiral expansion of the integral contains terms of 
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order q'^~^. The order of the infrared singular part follows from the counting 
rules at tree level, because it is generated by the region of integration, where 
the momenta flowing through the propagators are of the same order as in the 
tree level graphs. The remainder of the integration region does not contain 
infrared singularities and may thus be expanded in an ordinary Taylor series. 
An evaluation of the integral at its threshold at — s+ — (mjv + Af^)^ nicely 
shows the two parts 

The infrared singular part / is proportional to A/^^"^, while the remainder R is 
proportional to m^~^ and does therefore not contain a singularity at M^^ = 0, 
irrespective of the value of d. Since the regular contribution stems from a 
region where the variable of integration k is of order 0(111, f^), it violates the 
tree-level counting rules. 



3.2 Relation to HBCHFT 

The infrared singular part of the loop integrals can be obtained exchanging 
the order of the chiral expansion and the loop integration. To perform the 
expansion, both the residual momentum = — m^Vfj, and the loop mo- 
mentum k are counted as 0{q). Each term in the expansion of the integrand 
has a trivial dependence on the nucleon mass 111^ and does therefore not con- 
tain a regular part. Performing the loop integration one obtains the infrared 
part of the integral by summing the individual terms of the series. 

Let us illustrate this again for the self energy integral. In this case the 
expansion amounts to replacing 

1 1 (k- r)2 



TO^ — (P — k)^ 2mjv V ■ (k — r) {2mjv v ■ (k — r)}^ 

and the relativistic self-energy integral gets replaced by a sum of integrals 
which arise in HBCHPT (see fig. O): 



1 (fc~r)^("^i) 
k A/2 _ ^2 [2r7i„(t) - k-v r)]" 



Let us evaluate the integrals at = s+ = [M-ir+rar^Y by setting — Mt^v^. 
After expressing {k — r)^ = {k — M^v)'^ in terms of the denominators of the 
integrals and by using the fact that closed loops of static propagators vanish. 
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Figure 3. Internal line insertions. The double line denotes the heavy baryon propagator 
[2miv V ■ {k — r) — ie]~^ , the cross an insertion of (k — r)'^. 



we obtain the relation 

M 

US + )^ -In-l{s + ). (5) 

Since we have expanded the nucleon propagator, the integrals /„ depend only 
trivially on mw (they don't have a regular part) and dimensional analysis 
gives Ji = c(d)M^~^/m„. Therefore 

n—l n—1 

An explicit calculation of /i(s-)_) shows that the constant c(d) is indeed the 
same as in (||). 

In the case of thej^elf-energy integral it is possible to sum the kinematic 
insertions explicitly cM, but for a general integral, the above procedure would 
be rather clumsy. It however illustrates the relationship between HBCHPT 
and infrared regularization: The HBCHPT representation of the one nucleon 
amplitudes is obtained by first performing the chiral expansion and then the 
loop integration. The infrared regularized amplitudes can be obtained from 
HBCHPT by resumming the graphs with kinematic insertions. 

3.3 Feynman parameterization 

The explicit expression for the loop integral H{P^) in d dimensions is lengthy, 
but the splitting into the infrared singular part / and regular part R is easily 
obtained in the Schwinger-Feynman representation of the integral. 

, d'^k 1 1 f\ f d'^k 1 
H — — r — — — I dz 



{2nY a b Jo J (27r)'^ [ (1 - z) a + z 6 
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dz / 7 -7 7- ^ r-?^ = I + R . 6) 

J {2TTy [{l~z)a + zb]^ ^ ' 

The case of a general one loop integral (see fig.§b) can be reduced to the 
self-energy integral, by first using the ordinary Feynman parameterization to 
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combine all the meson propagators into a single one. After combining also 
the nucleon propagators in the same way, the general integral can be written 
as an integral over the self-energy integral H and derivatives with respect to 
its masses. Its infrared part is then obtained by replacing H by the Feynman 
parameterization for / given in . Feynman parameterizations of the infrared 
singular and regular parts of a general one loop integral can be found in the 
literature u. 

3.4 Chiral symmetry 

In the framework of the relativistic effective theory, the evaluation of an ampli- 
tude to one loop yields three contributions, arising from (a) the tree graphs, 

(b) the infrared singular and (c) the infrared regular part of the one loop 
graphs. At a given order in the chiral expansion, the one particle irreducible 
components of the regular part are polynomial in the external momenta. In 
coordinate space this contribution is given by local terms and thus equivalent 
to the tree graphs generated by a suitable Lagrangian AC. So that, if we 
replace the effective Lagrangian by 

£off = Ccs + AC , 

we may drop the regular part of the loop integrals. 

We show now that the modified Lagrangian C'^g occurring in our frame- 
work respects the Ward identities of chiral symmetry. In the limit of vanishing 
quark masses these identities interrelate the Green's functions of the form 

< N{P', S')\T {Oi{xi) . . . OniXn)} |iV(P, S) > , 

where the operators Oi stand for vector, axial- vector, scalar and pseudo-scalar 
quark-currents . 

Since Ccff is chirally symmetric, the tree level amplitudes (a) respect the 
Ward identities. Because dimensional regularization preserves the symmetry 
of >Cofr, the same is true for the loop contribution (b)+(c) for arbitrary values 
of the regularization parameter d. Since (b) involves only fractional and (c) 
only integer powers of the chiral expansion parameter, the sum of the two can 
only fulfill the Ward identities if they are separately fulfilled by the two parts 
and we conclude that the Lagrangian A>C, which generates the contributions 

(c) is indeed chirally symmetric. 

4 Representation of the TrA^-scattering amplitude at 0{q'^) 

The virtues of CHPT are twofold: One one hand chiral symmetry interrelates 
different amplitudes, on the other hand the framework allows one to obtain a 
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Figure 4. Kinematics of the elastic ttN scattering amplitude P, q {P',q') denote the mo- 
mentum of the incoming (outgoing) nuclcons and pions, and a (a') stands for the isospin 
indices of the incoming (outgoing) pion. 



parameterization of a given amplitude in a number of low energy parameters 
not fixed by the symmetry. 

In the following we focus on the energy dependence of the elastic pion 
nucleon scattering amplitude in the isospin lirtut, obtained on the basis of a 
calculation in infrared regularization to 0{q'^)n We explain, how one arrives 
at a representation of the amplitude in terms of functions of a single variable 
without invoking an expansion of its infrared singularities. As discussed in 
section ^, the chiral expansion of the triangle graph does not converge around 
its threshold. The expansion of the other graphs is convergent, but the con- 
vergence is in general rather slow and depends on the choice of the kinematic 
variables, which are kept fixed while one performs the expansion. 

In section |[ we have listed some of the analytic properties of the scattering 
amplitude. In the chiral representation of the amplitude, the corresponding 
cuts are described by loop integrals. These loop integrals in general also con- 
tain singularities outside the low energy region and we will in the following 
simplify the representation of the amplitude by chirally expanding the contri- 
butions from these singularities. As a result of these simplifications, we can 
represent the amplitude as the sum of the Born term, a polynomial contain- 
ing the low energy constants and nine functions of a single variable, which are 
given in terms of integrals over the imaginary parts of the loop integrals. 

In the context of our low energy analysis, the standard decomposition of 
the scattering amplitude into invariant amplitudes A and B is not well suited, 
because their leading order contributions cancel. Replacing Ahy D = A + i/B, 
with V — (s — u)/4mjv, the scattering amplitude is given by 
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where the upper index refers to the isospin decomposition 

The evaluation of the chiral perturbation series to one loop allows us to cal- 
culate the amplitudes and to 0{q^) and 0{q^), respectively. 

4-1 Simplification of the t-dependence of the amplitude 

The triangle graph with two nucleon propagators, fig. |l]d, and the box graph, 
fig. |l]e, have a cut for t > 4to^. Because this cut lies outside the low energy 
region, we can expand the corresponding graphs in the variable t. The higher 
order contributions in t are suppressed as t/Amj^ and are therefore also of 
higher order in the chiral expansion. After this simplification, the amplitude 
can be expressed in terms of functions of a single variable: The contribution of 
the graphs with an s-channel cut is a linear function of t and the contribution 
of the graphs with a t-channel cut is linear in v. 

4-2 Dispersive representation 

By construction, the infrared regular part of the loop integrals has no singular- 
ities in the low energypjegion. It however has unphysical singularities outside 
the low energy regionB and since the infrared part is given by the difference 
between the full integrals and their regular part, these singularities are also 
present in the infrared part. The infrared part of the self energy integral, for 
example, has a cut for s < and a pole at s = 0. Also these singularities can 
safely be chirally expanded. Let us illustrate in the example of the forward 
amplitude D^{v,t = 0) how this can be achieved (we use the bar to indicate 
that the pseudo- vector Born term has been subtracted from D^). In a first 
step, we expand the amplitude around ly — to 0(z^*): 

D+{iy,0)=do + dy- + d2iy^ + 0{iy^) 

The absence of odd powers of v is due to crossing symmetry. Note that we 
have not expanded the coefficients di in powers of the quark masses, because 
this expansion converges only slowly, especially in the case of c?2- Once its con- 
tributions to the polynomial part are subtracted, the remaining contribution 
from the physical s- and u-channel cut is given by the dispersion integral 

2 r°° dv' 
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lmsD^{i', 0) denotes the discontinuity over the right hand cut of the s-channel 
graphs, which is a consequence of unitarity and is the same in dimensional 
and infrared regularization. 

If we subtract from the amphtude both the polynomial and the dispersion 
integral, we are left with 

D+^{i^) = D+{i^, 0) - D+{:^) - do - diiy^ - d2i^^ , (8) 

which receives contributions only from the singularities outside the low en- 
ergy region. Since we have accounted for all infrared singularities and since 
-^reg(^jO) proportional to v^, its chiral expansion involves only terms of 
order 0{q^) and higher and we can drop D^^^ to the accuracy of our calcula- 
tion. The representation of the forward amplitude now takes the form 

D+{,y,0)^do + diiy^ +d2iy'^ + D+{i^) + 0{p^). (9) 

The coefficients di contain low energy constants, whose values are not fixed by 
the symmetry. By construction D^{v) does not contribute to the coefficients 
di, so that we can identify these constants with tha^lowest coefficients of the 
subthreshold-expansion of Hohler and coUaboratorsH. 

We can follow the same strategy to obtain a dispersive representation 
of for nonzero values of t. This representation involves three more sub- 
traction constants, corresponding to the coefficients of <, iP t and of the 
subthreshold expansion, as well as a dispersion integral over the t-channel 



cut. After the simplifications made in section 4.1, the s-channel imaginary 
part is linear in t and one needs an additional dispersion integral, which sup- 
plies the contributions of the s- and u-channel cut which are proportional to t. 
In the case of the amplitudes -B^, the s-channel imaginary part is independent 
of t, so that one function suffices to describe the s-channel cut contribution 
for each of these. Furthermore the i-channel imaginary part of _B+ vanishes 
to the order of our calculation, so that we need a total of nine functions of a 
single variable to account for all the low energy singularities of the amplitude 
to 0{q^). 



5 Conclusion 



We have presented a formulation for BCHPT (infrared regularization) in 
which both, Lorentz invariance and the chiral power counting rules are pre- 
served at every stage of the calculation. The method relies on the fact that 
for noninteger values of the regularization parameter d, the infrared singular 
part of the relativistic loop integrals can be unambiguously separated from 
the remainder. Since the remainder is free of infrared singularities, it is a 
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polynomial in the external momenta and quark masses to any finite order of 
the chiral expansion and can be absorbed into a redefinition of the coupling 
constants of the effective Lagrangian. 

We have shown that the chiral expansion of the tt TV-scattering amplitude 
does not converge in the vicinity of the threshold at t = 4M^. The breakdown 
of the expansion is related to the fact that the expansion of the nucleon 
propagator around the static limit destroys some of the analytic properties 
of the amplitude. In the relativistic formulation of the effective theory the 
correct analytic structure is guaranteed and one obtains a representation, 
which is valid in the entire low energy region. 

The representation of a given amplitude obtained in the relativistic frame- 
work contains in general also singularities outside the low energy region. We 
have shown how to simplify the representation of the TriV-scattering ampli- 
tude at 0{q'^), by chirally expanding the contributions of these singularities. 
As a result of these simplifications, the scattering amplitude can be written 
as the Born term, a polynomial containing the low energy constants and nine 
functions of a single variable. The corresponding functions are given in terms 
of dispersion integrals over the imaginary parts of the one loop graphs. 



References 



1. E. Jenkins and A. V. Manohar, Phys. Lett. B 255, 558 (1991). 

V. Bernard, N. Kaiser, J. Kambor and U.-G. Meissner, Nud. Phys. B 
388, 315 (1992). 

2. J. Gasser, M. E. Sainio and A. Svarc, Nud. Phys. B 307, 779 (1988). 

3. T. Becher and H. Leutwyler, Eur. Phys. J. C 9, 643 (1999) 

4. H.-B. Tang, A new approach to chiral perturbation theory for matter 
fields, |hcp-ph/9607436| 

P. J. Ellis and H.-B. Tang, Phys. Rev. C 57, 3356 (1998) 

5. J. Schweizer, Low energy representation for the axial form factor of the 
nucleon, diploma thesis, Bern 2000. 

6. B. Kubis and U. Meis sner, Low energy analysis of the nucleon electro- 
magnetic form factors, hcp-ph /0007056 . 

B. Kubis and U. Meissner, Baryon form factors in chiral perturbation 



theory, |hep-ph/0010283 . 

T. Becher and H. Leutwyler, to appear. 

G. Hohler, in Landolt-Bornstein, 9b2, ed. H. Schopper (Springer, Berlin, 
1983). 



12 



